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2^ ■ Abstract 

^ . We investigate the spin-magnetism of mesoscopic metallic grains. In the av- 

^^ , erage response of an ensemble of grains there are corrections to macroscopic 

behaviour due to both spectral fluctuations and electron-electron interactions. 

These corrections are a non-linear function of the magnetic field. Their tem- 

perature dependence is calculated numerically and analytically. An experi- 
(^ . ment is proposed to measure the unknown interaction coupling constant in the 

CS , cooper channel. For a single sample the magnetization is found to fluctuate 

reproducibly about the mean. These fluctuations directly probe the energy 

level statistics. 
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In the early days of quantum statistical mechanics it was realized by Frohlich that the 
thermodynamic properties of an ideal fermi gas in a finite box would change qualitatively 
from its macroscopic behavior for temperatures (kT) much smaller than the mean level 
spacing (A^) [l|]. Frohlich predicted that the large level-spacing (compared to kT) would 
lead to exponential temperature dependence of thermodynamic quantities ~ exp(— Ae/fcT) 
in normal metals in contrast to the weak T-dependence found in bulk metals. Subsequently 
Kubo and Gorkov and Eliashberg |^ took into account the fluctuations in the level 
spacing of metallic grains and found that upon averaging over an ensemble of grains the 
exponential dependences were softened to power laws which depended on the precise level 
statistics. Considerable effort has been expended in experimental searches for such finite 
size corrections to thermodynamics in clusters of small metallic particles @] using NMR and 
other techniques. 

Recent work in mesoscopic physics has pointed out that there is another relevant energy 
scale for the thermodynamics of small particles, the Thouless energy E^ = h/ru, the inverse 
time to diffuse across the sample. In a micron-size metal film Ef. ~ IK ~ lO^A^ so it is 
straightforward to perform experiments under the conditions Ec ^ kT ^ A^, which we 
will term the mesoscopic regime. Work on orbital magnetism has shown that in the meso- 
scopic regime corrections to the thermodynamics are determined by the long-range spectral 
fluctuations, not the level spacing distribution. However the recent work on mesoscopic 
thermodynamics has tended to focus on the orbital magnetization and the associated per- 
sistent currents and has not made contact with the work cited above which focused on spin 
magnetization (and speciflc heat). That early work ignored the effect of long-range spectral 
fluctuations and found exponential convergence to macroscopic behavior for kT 3> A^ |p. 

In this Letter we determine the effect of spectral fluctuations on the spin magnetization. 
The spin magnetization of course contributes to the experimentally measured magnetization 
and must be understood if it is to be separated from the orbital magnetization. Its quantum 
corrections have been well-studied in the microscopic {kT <^ A^) and macroscopic regimes 
0J^. We flnd that in the mesoscopic regime the quantum corrections to the average spin 
magnetization are not exponentially small in kT/A^ and can be comparable to the correc- 
tions to the orbital magnetization. There are corrections due both to interaction effects and 
to spectral fluctuations. Experiments that probe the thermodynamics of mesoscopic sys- 
tems have been performed on both individual samples [^ and arrays and thus both the 
mean and variance of the magnetization are in principle measurable. Our theory indicates 
that the interaction corrections will give the dominant contribution to the mean while the 
long-range spectral fluctuations will be measurable in the variance, providing a direct probe 
of Wigner-Dyson level statistics |^. The interaction corrections to the mean will allow the 
determination of a crucial interaction constant which also controls the size of the persistent 
current. 

Initially we neglect electron-electron interactions and calculate the corrections to the 
Pauli susceptibility due to spectral fluctuations in the mesoscopic regime. In typical exper- 
iments, where the sample dimensions are of order a micron or smaller, the thermodynamic 
properties should correspond to the canonical ensemble (CE) since the electron number N 
on each specimen is flxed by charge neutrality 0. A standard procedure in statistical me- 
chanics is to replace a canonical average with the average in the equivalent grand canonical 
ensemble (EGCE), which is a grand canonical ensemble (GCE) with the chemical potential 



adjusted to have an average number equal to N. Imry pointed out that this adjustment led 
to new quantum interference contributions to thermodynamics proportional to the density 
of states (DOS) fluctuations |T^. These average to zero in macroscopic systems but become 



important in mesoscopic systems; hence Imry's approach has become the standard technique 
for calculating mesoscopic corrections to thermodynamics ^^. We note however that the 



CE and the EGCE are only equivalent in the limit A^ ^ 1 and kT ^ A^ ||12|- However 
in this work we find by direct numerical simulations that the finite-size deviations between 
the CE and EGCE are negligible in the mesoscopic regime; so we will proceed using Imry's 
method. In this method the leading correction to the free energy of the EGCE is given by 

SFm = ^ dE dE' < 5p{E)5p{E') > f{E)f{E'). (1) 
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Here f{E) denotes the fermi function, Sp{E) is the fluctuating part of the DOS, and < . . . > 
denotes an average over impurity realizations. 

In order to compute the magnetization corrections using Eq. (1) we must obtain the 
magnetic field dependence of 6Fi\i. Here we focus on systems with negligible spin-orbit 
interaction; in this case the DOS fluctuation may written as 

6po{E + A,/2) + 6po{E - A,/2). (2) 

Here A^ denotes the Zeeman splitting and 6pq{E) denotes the DOS fluctuations for spinless 
electrons in a magnetic field. If this quantity is known the dependence of SF^ on A^ and 
hence the correction to the spin magnetization follows by inserting Eq. (2) into Eq. (1) and 
differentiating. 

It is now established by microscopic calculations that the energy levels of disordered 
metallic samples are well-described by Wigner-Dyson random matrix theory(RMT) [Q for 



energy intervals up to Ec U3|,|14[. We find that the deviation from RMT does not affect 
the spin magnetization when kT -C Ec and hence we may assume that the DOS correla- 
tion function in Eq. (1) is that given by RMT. This correlation function decays slowly as 
1/i3[7t{E — E')]"^, for \E — E'\ much larger than A^ [|l^,|l^. Here (3 = 1,2 depending on 
whether or not the orbital effects of the magnetic field are sufficient to break time reversal 
symmetry. Using this form in Eq. (1) and differentiating one finds 

AM = -,.-5|(2rtr)i:Re|^— ^} (3) 

where pe denotes the electron's magnetic moment, z/ is a bosonic Matsubara frequency, and 
the energy integrals have been converted to Matsubara sums. This result is insensitive to 
the cut-off of the sums in the regime of interest, Ec ^ kT ^ A^. The sum as a function of 
A^ is easily evaluated numerically, and analytically in the limits of small and large Zeeman 
splitting. 

AM = ^C(3)/ie(A,A,)/(A;T)2 A, « kT; 

= 2^/ie(A,/A,) A, > kT. (4) 



At low field the correction to the magnetization is linear and decreases as (Ag/fcT)^; thus 
as expected it is not exponentially small as found when long-range energy level fluctuations 
are completely neglected. Furthermore, the magnetization is non-linear at high fields, for 
Zeeman splittings of the order oikT. The Pauli magnetization, Mp = 2fieAz/A^, by contrast 
remains strictly linear at all relevant fields. Thus despite its small size compared to Mp, the 
correction may be distinguished by its non-linear susceptibility ||^. At T ~ 0.1 K the field 
scale at which the non-linearity sets in is i? ~ 0.1 T, so it occurs at experimentally accessible 
parameter values. To emphasize that the (non-interacting) spin magnetization correction 
probes the long-range level fluctuations one can easily evaluate it for the (Poisson) case 
of uncorrelated levels. Since the fluctuations are enhanced one finds a strongly enhanced 
AM ~ {A^/kT)Mp. It is possible that this behavior may be observable in appropriate 
ballistic semiconductor microstructures, but this question requires further study. 

We now discuss the variance of the magnetization (which arises due to quantum mechani- 
cal coherence). Following Imry's method it can be shown that chemical potential adjustment 
does not affect the variance to second order in 5p, so one may employ the more convenient 
GCE. Using Eq. (2) the variance of the magnetization is given by 
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5M' =< I J dE Y: a6po{E + ^)f{E)] >, (5) 

so again the mesoscopic corrections are determined by the correlation function of the DOS. 
Using the RMT correlation function in Eq. (5) leads to 

6M' = ^C(3)/i^(A,/fcT)2 ^^ ^ ^^. 

1 pl\n{A,/kT) A, >A;T. (6) 
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Note that the root-mean-squared magnetization fluctuation is parametrically larger than the 
correction to the mean. It grows linearly with Zeeman splitting at low field and slows to a 
logarithmic increase at high fields. The behavior for kT ^ A^ follows from our previous re- 
sult for AM within Imry's method; if one expands Eq. (5) to 2nd order in A^ and integrates 
by parts one finds AM/Mp = 6M'^/Mp, independent of the nature of the level statistics. It 
thus follows that in general 5M/AM = (Mp/AM)^/^ 3> 1, and the magnetization fluctua- 
tions at low field will always be much larger than the correction to mean. The saturation of 
5M^ to a logarithmic increase at high-fields (A^ ^ kT) also has a simple interpretation. At 
T=0 the GCE magnetization literally counts the number of levels (at zero splitting) in an 
energy window of width A^ centred at e^. Thus 6M'^ measures the variance of that number. 
It is a famous result of RMT that this variance increases logarithmically with the size of the 
interval Q] (as opposed to the linear increase for uncorrelated levels); hence the logarithmic 
dependence on A^. 

By the familiar "ergodic hypothesis" of mesoscopic physics |T^ we expect the statistical 



fluctuations of M to manifest themselves as fluctuations of the magnetization of a given 
specimen as a function of external parameters (e.g. fermi energy or magnetic flux). However 
simply varying magnetic field is not a good method as this changes A^ as well. With a 2D 
mesoscopic sample one can envision varying the tilt angle of the field with respect to the 
plane of the sample, hence varying the magnetic flux at fixed A^. The main difficulty with 
such an experiment is the small absolute size of the effect, which is of order a Bohr magneton 



per sample. However NMR has the sensitivity to measure such small effects (roughly 1-10 
ppm for a micron-size film at i? ~ IT) if appropriate samples can be obtained. Here the 
mesoscopic size scale offers entirely new possibilities as one can imagine fabricating arrays of 
metal samples lithographically (as was done to measure persistent currents [§]) and achieving 
much greater homogeneity than in conventional NMR of small particles. The mesoscopic 
fluctuations would appear as a broadening of the NMR line which increases with decreasing 
temperature 0. 

As noted above, Imry's method assumes that the EGCE, which uses fermi functions 
with an adjusted chemical potential yields, the same thermodynamic properties as the true 
canonical ensemble (which can never be exactly represented by a fermi function), at least 
in the mesoscopic regime {Ec ^ kT ^ A^). Since it has been shown that the CE and 
EGCE have substantially different thermodynamic properties when kT < A^ [Q, it seemed 
necessary to test their equivalence in the mesoscopic regime. Moreover Imry's formula 
(Eq. (1)) is the first term in an expansion involving successively higher level correlation 
functions; therefore its quantitative accuracy for a particular range of values of kT/A^ needs 
to be tested. 

To do this we have calculated the exact ensemble-averaged CE and EGCE spin magneti- 
zation numerically. The ensemble of level-sequences was generated by direct diagonalization 
of 10^ 1000 X 1000 gaussian orthogonal (/3 = 1) random matrices and appropriate "unfold- 
ing" of the levels obtained @jl3l- Evaluation of the CE partition function (for a given level 



sequence) by brute force summation is not possible for kT ^ A^ due to the enormous num- 
ber of relevant states. To bypass these difficulties we use the Darwin- Fowler representation 
of the CE partition function [0]. In this approach an infinite (Darwin- Fowler) polynomial 
is defined whose coefficients are precisely the partition function for different values of A^. 
Each coefficient can be extracted (in principle) by an appropriate contour integral. Eval- 



uation of this integral by steepest descent yields the EGCE partition function [T^. We 
have shown that it is possible to represent the magnetization directly as a Darwin-Fowler 
type contour integral, which can be evaluated numerically for sequences of 500 levels with- 
out difficulty. Exact numerical calculations within the EGCE are straightforward as one 
may use the fer mi-function and simply adjust the chemical potential appropriately for each 
level-sequence. This procedure includes the higher order correlations neglected in Eq. (1). 

The results of these calculations for the magnetization are shown in Fig. 1. They may 
be summarized as follows: 1) When kT > A^ the CE and EGCE averages agree almost 
perfectly. Hence the EGCE is an excellent approximation in the mesoscopic regime [IT^ 



2) The agreement between the numerical and analytic calculations of SM"^ is very good 
whether the variance is calculated in the (unadjusted) GCE or the EGCE. 3) The analytic 
formula for the mean correction AM and the numerical calculations differ substantially in 
the interval A^ < kT < 10 A^ (roughly by a factor three). It is also not possible to fit the 
numerical results to the quadratic T-dependence predicted (although this may be due to the 
breakdown of the accuracy of our numerics at the highest values of kT). This shows that 
there are substantial corrections to Eq. (1) in the range of parameters where the mesoscopic 
corrections are largest ( A^ < kT < 10 AJ. 

Finally we include the effects of interaction on the mean spin magnetization using per- 
turbation theory. It is known that these effects are important for persistent currents in 
disordered metals. The relevant diagrams for the spin magnetization are already known in 



the literature pO[, only their behavior for kT < E^ had not been explored. We find 



AM.„.Hf + A.)*ri:i:Re|j^;^^^=f^^| (7) 

where F, Ac are coupling constants in the diffuson, cooperon channels which we will regard 
as parameters to be determined experimentally. Here Vt^ denotes the eigenvalues of L^V^ 
(L is a typical sample dimension) with Neumann boundary conditions. ^2^ = is always an 
eigenvalue and it is easy to show that this term dominates the sum when kT < E^. Keeping 
only this term and performing the Matsubara sums yields 

AM,^t = {F + X,)j,a2){AJkT) A, «A;T; 

F + Xc A^ > kT. (8) 

This correction is parametrically larger than that predicted by the non-interacting theory 
(Eq. 4) if the dimensionless interaction constants are not too small. Note that AMi^t 
is linear at low field (A^ -C kT) and saturates at high field (A^ ^ kT). It is therefore 
highly non-linear for fields around A^ ~ kT. In this respect its response mimics that of 
an isolated spin which may complicate its experimental determination. However because 
it involves phase-coherent diffusion the interaction correction is very anisotropic. Eq. (8) 
assumes that the cooperon term is not suppressed; the magnetic fiux at which this term 
is suppressed is the fiux quantum h/e which corresponds to a normal magnetic field of 
order 10^ G. As noted above, if samples are fabricated as small metal films the fiux may be 
varied by tilting the sample allowing a measurement of the term proportional to Ac- This 
would be an important result of such an experiment as Ac controls the size of the persistent 
current and its value is only measurable very indirectly in non-superconducting materials 
PT| . If such a measurement, e.g. in copper, agreed with the measured amplitude of the 



persistent current [g| it would provide strong evidence that the interaction correction is the 
relevant one in that case as well. Finally, in the case of persistent currents the interacting 
and non-interacting theories both give exponential temperature-dependence whereas for the 
spin magnetization different power-laws are predicted (see Eqs. (4), (8)) which may make it 
easier to experimentally distinguish the two effects. 

To summarize, we have studied the field and temperature dependences of mesoscopic 
corrections to macroscopic spin-magnetism for grains of normal metal. Measurement of the 
correction for an ensemble of grains would provide an independent estimate of an important 
interaction coupling constant which is difficult to measure in non-superconducting metals. 
Single samples are found to exhibit magnetization fiuctuations that depend only on the 
external parameters, the applied field and temperature, not on sample properties such as 
size, impurity density or fermi energy. 

We thank H. Baranger and S. Barrett for helpful comments. This work was partially 
supported by NSF grant DMR-92 15065. 
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FIGURES 

Comparison of the temperature dependence of the average magnetization at low field, 
calculated in the different ensembles: Squares represent the numerical CE result; triangles, 
the numerical EGCE result; and the smooth curve, the analytical result of Eq (4). Inset: 
The same results on a log-log plot. The magnetization fluctuations (normalized by Mp) are 
also shown. With this normalization, the analytical curve for the fluctuations coincides with 
that for the average (see text); circles represent the numerical GCE result. 



